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LINEAR  THEORY  OF  A  QUASIOPTICAL  GYROKLYSTRON 
WITH  NONUNIFORM  MAGNETIC  FIELD 


I.  Introduction 

Gyrotrons  are  currently  under  development  as  high  average-power  sources  of  millimeter  waves 
for  electron  cyclotron  resonance  heating  (ECRH)  of  fusion  plasmas.  The  leading  candidate  for 
ECRH  is  the  waveguide  cavity  gyrotron,  which  has  achieved  powers  of  1.2  MW  with  20%  efficiency 
at  a  frequency  of  148  GHz  in  microsecond-length  pulses.1  An  alternative  configuration  is  the 
quasioptical  gyrotron  (QOG),2  which  has  several  advantages  over  the  conventional  cavity  gyrotron 
as  the  power  and  frequency  of  the  radiation  are  increased. 

The  resonantor  of  the  QOG  is  formed  by  a  pair  of  spherical  mirrors  separated  by  many 
radiation  wavelengths.  The  device  can  be  made  to  oscillate  in  the  fundamental  transverse  mode 
only,  since  higher-order  modes  suffer  much  larger  diffraction  losses  due  to  the  finite  size  of  the 
mirrors.  The  axis  of  the  open  resonator  is  oriented  perpendicular  to  the  path  of  the  electron  beam, 
so  that  the  radiation  is  collected  separately  from  the  beam.  This  separation  is  advantageous  for 
high  power  operation  and  facilitates  the  use  of  depressed  collector  technology  for  higher  device 
efficiency.  The  peak  ohmic  heating  density  on  the  mirrors,  which  is  important  for  continuous- wave 
(cw)  devices,  is  reduced  to  an  allowed  level  by  increasing  the  separation  between  the  resonator 
mirrors.  Increasing  the  mirror  separation  for  a  given  set  of  mirrors  increases  the  output  coupling 
of  the  radiation,  although  the  radiation  waist  at  the  center  of  the  resonator  is  insensitive  to 
changes  in  resonator  length.  In  this  manner  the  output  coupling  can  be  varied  while  the  beam- 
wave  interaction  length  remains  approximately  constant.  This  feature  is  unique  to  the  QOG  and 
provides  a  technique  to  optimize  the  strength  of  the  electric  field  in  the  resonator. 

The  axial  mode  separation  in  the  QOG  is  small  compared  to  the  interaction  bandwidth  in  cw- 
relevant  devices,  therefore  multimode  effects  are  important.  Most  gyrotron  applications  require 
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a  single-mode  output,  so  that  suppressing  unwanted  longitudinal  modes  is  of  interest  in  the 
QOG.  Even  if  multimode  operation  is  permitted,  it  is  not  dear  whether  the  optimum  efficiency 
of  single-mode  operation  can  be  obtained.  Multimode  theoretical  effidendes  are  comparable  to 
optimum  single-mode  effidendes  in  the  QOG,3,4  although  both  of  these  predictions  overestimate 
experimental  results  to  date. 

Recent  theoretical  work  has  shown  that  there  exist  large  regions  in  operating  parameter  space 
where  the  output  is  single-moded  even  though  the  axial  mode  density  is  high.5  Single-mode 
operation  has  been  observed  up  to  250  kW  in  a  cw-relevant  QOG  at  120  GHz  with  a  13  psec 
pulse  length.6  However,  the  gyrotron  had  to  be  detuned  slightly  to  obtain  single-mode  operation, 
which  reduced  the  output  efficiency  from  the  multimode  value.  The  output  was  multimoded  at 
higher  powers  with  as  many  as  5  or  6  longitudinal  modes  present  at  500  kW  output  power.  The 
region  in  parameter  space  of  single-mode  operation  in  the  experiement  was  much  less  than  that 
predicted  by  the  theory.  Theoretical  predictions  also  indicate  that  tilting  the  axis  of  the  resonator 
a  few  degrees  should  increase  the  regions  of  stable  single-mode  operation  for  an  annular  dectron 
beam.  Experimental  results  indicate  that  tilting  the  resonator  improves  the  single-mode  efficiency, 
although  the  region  of  single-mode  operation  does  not  increase  significantly.  6 

One  technique  which  can  be  used  to  select  the  optimum  longitudinal  mode  in  the  output 
resonator  is  to  use  an  external  source  to  set  up  the  rf  field  at  the  desired  frequency.  This  may 
be  accomplished  by  direct  injection  of  radiation  into  the  output  resonator7,8  or  by  prebunching 
the  electron  beam  in  an  upstream  resonator.8  Direct  injection  requires  a  high  power  circulator 
or  isolator  so  that  the  external  source  is  not  locked  by  the  gyrotron.  Such  an  isolator  is  difficult 
to  obtain  at  millimeter  wavelengths,  so  prebunching  the  electron  beam  is  considered  here.  A 
two-resonator  quasioptical  gyroklystron,  shown  in  Figure  1,  comprises  a  pair  of  open  resonators 
separated  by  a  field-free  drift  region.  A  low  amplitude  electric  field  is  excited  in  the  first  resonator 
via  an  external  millimeter-wave  source,  such  as  an  extended  interaction  oscillator  (EIO).  Electron 
cyclotron  absorption  provides  electrons  in  the  prebunching  resonator  with  a  “kick”  in  perpendicu¬ 
lar  momentum,  depending  upon  their  entrance  gyrophase,  so  that  the  particles  gyrate  at  different 
frequencies  in  the  drift  region.  If  the  strength  of  the  electric  field  in  the  first  resonator  and  the 
length  of  the  drift  region  are  chosen  properly,  the  electrons  will  arrive  at  the  output  resonator 
strongly  bunched  in  gyrophase  angle  and  will  efficiently  give  up  their  energy  to  the  wave  fields. 

In  this  work  we  consider  the  case  where  the  beam  current  is  above  the  threshhold  current  of 
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the  output  resonator  in  the  absence  of  a  prebunching  signal.  Thus,  this  regime  of  operation  is 
fundamentally  different  than  the  gyroklystron  amplifier.  The  strength  of  the  prebunching  signal 
and  the  frequency  difference  between  the  two  sources  has  a  profound  effect  on  the  operation  of  the 
gyroklystron.  If  the  prebunching  signal  is  sufficiently  large  and  the  frequency  separation  small,  the 
gyroklystron  will  operate  as  a  locked  oscillator  where  the  equilibrium  has  a  specified  phase.9  Here 
the  steady  state  operation  of  the  output  resonator  depends  upon  the  amplitude,  frequency,  and 
phase  of  the  prebunching  signal.  It  is  in  this  regime  that  the  electrons  are  strongly  prebunched 
and  the  peak  gyroklystron  transverse  efficiency  is  higher  than  that  of  the  single  cavity  gyrotron.10 
The  drive  signal  is  present  throughout  the  duration  of  the  oscillation  and  has  a  strong  effect  on 
the  steady  state  of  the  oscillator. 

A  second  region  of  operation  is  the  mode-primed  gyroklystron,  where  the  external  source 
is  used  only  during  the  build-up  of  fields  in  the  output  resonator.8  The  prebunched  electron 
beam  selects  a  particular  longitudinal  mode  in  the  output  resonator,  which  is  given  an  initial 
advantage  over  neighboring  modes.  If  this  mode  is  stable  with  respect  to  decay  into  sidebands, 
it  will  grow  and  nonlinearly  suppress  satellite  modes.  The  prebunching  signal  is  used  to  mode- 
prime  the  oscillator,  and  does  not  affect  the  frequency  of  the  final  state.  The  main  advantage  to 
mode-priming  is  that  the  beam  premodulation  need  only  be  large  compared  to  the  noise  present 
at  the  start  of  oscillation  in  the  output  resonator.  This  reduces  the  power  requirement  for  the 
millimeter-wave  source,  which  is  often  a  limitation  at  high  frequencies. 

Another  configuration  is  the  mode-locked  gyroklystron,  where  a  portion  of  the  output  is  fed 
back  into  the  input  resonator.11  The  mirror  separation  of  the  prebunching  resonator  is  chosen 
so  that  only  one  mode  in  the  output  spectrum  is  resonant  in  the  prebuncher.  Thus,  the  input 
resonator  acts  like  a  filter  and  the  feedback  of  the  desired  longitudinal  mode  is  accomplished 
without  an  external  source.  In  the  mode-locked  case  the  frequency  of  the  oscillator  is  fixed  and 
there  is  a  relative  phase  between  the  two  resonators,  although  the  equilibrium  doesn’t  have  a 
selected  phase. 

An  alternative  technique  which  does  not  require  an  external  source  is  the  gyroklystron  oscillator 
This  device  is  operated  at  a  beam  current  slightly  above  the  threshhold  current  for  the  prebunch¬ 
ing  resonator.  The  electric  field  in  the  first  resonator  will  modulate  the  electron  beam  and  excite 
that  frequency  in  the  output  resonator.  If  the  current  is  increased  further,  the  strong  rf  fields  in 
the  input  resonator  will  introduce  a  large  energy  spread  to  the  electrons  and  spoil  the  efficient 
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oscillation  in  the  output  resonator.  Thus,  this  type  of  device  wall  only  operate  over  a  limited 
range  of  parameters. 

Two  major  limitations  on  conventional  cavity  gyroklystron  experiments  are  spurious  oscilla¬ 
tions  in  the  input  cavity,  output  cavity,  and  drift  region  and  velocity  spread  of  the  electron  beam. 
The  QOG  has  a  distinct  advantage  with  regard  to  spurious  oscillations  because  modes  other  than 
the  TEMoov  are  strongly  discriminated  against  due  to  the  finite  size  of  the  mirrors.  Also,  only 
a  single  longitudinal  mode  is  within  the  gyrotron  interaction  bandwidth  for  small  mirror  separa¬ 
tions.  Hence,  exciting  the  proper  mode  in  the  input  resonator  should  be  straightforward.  Velocity 
spread  of  the  beam  is  caused  predominantly  by  the  electron  gun  and  is  difficult  to  measure  and 
control.  Thermal  spread  manifests  itself  by  debunching  the  electrons  as  they  drift  between  the 
two  resonators.  The  most  common  means  of  addressing  this  problem  is  to  minimize  the  length  of 
the  drift  region,  which  sets  a  limit  on  the  maximum  gyrophase  bunching. 

The  QOG  experiment  at  the  Naval  Research  Laboratory  uses  a  superconducting  magnet  with 
a  four-inch  crossbore  where  the  field  is  produced  by  a  pair  of  modified  Helmholtz  coils.  The  axial 
magnetic  field  has  a  7%  dip  at  the  center  of  the  crossbore,  which  is  the  location  of  the  resonator  in 
previous  experiments.  For  the  proposed  gyroklystron  experiment,  the  output  resonator  will  remain 
at  the  center  of  the  crossbore  while  the  prebunching  resonator  will  be  placed  upstream,  before 
the  magnetic  field  peak,  at  approximately  the  same  field  as  the  output  resonator.  The  frequency 
of  the  experiment  is  determined  by  the  availability  of  a  suitable  millimeter-wave  source,  in  this 
case  an  85  GHz  EIO.  The  EIO  generates  1.5  kW  peak  power  with  a  pulse  length  up  to  2  psec  and 
a  mechanical  tuning  range  of  ±1  GHz. 

This  paper  is  organized  as  follows.  Section  II  analyzes  the  effect  of  a  nonuniform  magnetic  field 
on  the  prebunching  of  the  electrons  in  the  first  resonator.  It  is  found  that  afield  gradient  reduces 
the  electron  bunching  parameter  somewhat  but  poses  no  serious  limitation  to  the  experiment. 
Section  III  examines  the  effect  of  velocity  spread  on  the  linear  efficiency  of  the  gy  roklystron.  The 
nonuniform  magnetic  field  can  be  used  to  minimize  the  deleterious  effect  of  thermal  spread  by 
detuning  the  cyclotron  frequencies  in  the  two  resonators.  This  calculation  determines  the  position 
of  the  prebunching  resonator  so  that  the  bunching  is  preserved  in  the  presence  of  velocity  spread  in 
the  drift  region.  Section  IV  provides  a  design  example  where  the  bunching  parameter  is  calculated 
for  the  experimental  parameters,  while  Section  V  contains  conclusions. 
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II.  Linear  Theory  in  the  Prebunching  Resonator 
A.  Uniform  Magnetic  Field 

We  first  consider  the  case  of  a  uniform  magnetic  field  and  calculate  the  bunching  of  the 
electrons  due  to  the  electric  field  in  the  input  resonator  and  ballistic  bunching  in  the  drift  region. 
The  equations  which  describe  the  spacial  evolution  of  the  slow  phase  and  the  perpendicular 
momentum  of  the  particles  can  be  written3 

+  (1) 

=  -^-E(z)smky,zos».  (2) 

These  equations,  known  as  the  slow  time  scale  equations  of  motion,  are  derived  under  the  as¬ 
sumption  that  the  wave  fields  evolve  on  a  time  scale  much  longer  than  the  cyclotron  period.  The 
beam-wave  interaction  is  taken  to  occur  at  the  fundamental  cyclotron  frequency  under  the  weakly 
relativistic  limit  (7  —  1  ■<  1).  The  slowly  varying  gyrophase  angle  is  defined  9  —  ip  -  u>t,  where 
is  the  particle  phase,  k  is  the  free  space  wavenumber,  yg  =  y  +  ( Px/raft )  cos  ip,  and  the  electric 
field  is  polarized  in  the  x-direction.  The  slow  phase  of  an  unperturbed  particle  is 

B  =  90  +  —(ft  -  7 w)*,  (3) 

Pz 

where  z  =  0  is  the  position  of  the  prebunching  resonator  and  9q  is  the  random  phase  of  the  particle 
at  the  entrance  of  the  resonator.  In  general,  the  relativistic  mass  factor  (7),  the  nonrelativistic 
cyclotron  frequency  (ft),  and  the  parallel  momentum  of  the  particle  (p2)  are  also  functions  of 
position. 

The  mirrors  which  make  up  the  resonator  are  chosen  to  support  only  the  fundamental  trans¬ 
verse  mode  (TEMoo).  Higher  order  transverse  modes  have  much  larger  diffraction  losses  and 
typically  have  quality  factors  three  to  five  times  less  than  the  fundamental  mode.  The  assumed 
radiation  profile  is  Gaussian 

E{z)  =  £01  exp  [-A-].  (4) 

Wqi 

Here,  uipi  is  the  radiation  waist  at  the  center  of  a  symmetric  resonator  and  the  subscript  denotes 
the  first  resonator.  The  radiation  waist  at  the  center  of  the  resonator  is  given  by13 


In  this  equation,  d  is  the  mirror  separation  and  Rc  is  the  radius  of  curvature  of  the  mirrors. 

5 


The  change  in  perpendicular  momentum  for  an  electron  which  traverses  the  prebunching 
resonator  with  a  low  amplitude  electric  field  can  be  calculated  by  integrating  Eq.  (2). 

Apj.  =  — f  E(z)8inkygcos[60  +  —(ft  -  70 u)z]dz.  (6) 

2p*o  J —00  Pzo 

Here  it  has  been  assumed  that  pz  =  pz 0  is  a  constant  of  the  motion  and  the  relativistic  mass  factor 
7  is  constant  through  the  first  resonator.  If  the  magnetic  field  is  uniform  through  the  prebunching 
resonator,  ft  =  fto  and  the  integral  can  be  performed  analytically  yielding 

\  y/nem-rowoi  n  .  ,  „  r-(ft0  -  7ow)2tn^m21 

Apx  = - — - -Eoismfcyflcos0oexp[— - jp - ].  (7) 

The  change  in  perpendicular  momentum  is  a  function  of  the  entrance  phase  of  the  particle,  so  that 
some  electrons  gain  energy  while  some  lose  energy.  Note  that  the  change  in  perpendicular  mo¬ 
mentum  is  due  to  the  electric  field  in  the  prebuncher  and  occurs  over  a  distance  of  approximately 

2u>oi. 

In  the  field-free  drift  region, 

f  =  (8) 

Here,  the  relativistic  mass  factor  has  been  modified  to  include  the  “kick”  from  the  prebunching 
re  "orator. 

Ta,*1+$$3>  (9' 

Thus,  the  slow  phase  of  the  particle  at  the  entrance  of  the  second  resonator  is  6  =  6q  +  A0  where 

m 

A 9  =  — (fto  -  7ow)(I  +  z)  4-  q  cos  0O  sin  k v-  (10) 

PzO 


y/KupxeEoirv01L  ,  (fto  -  7ow)2tt&m21 

,= - - expI - 3 S - 1'  (U) 

The  variable  z  has  been  redefined  so  that  z  —  0  is  now  at  the  center  of  the  power  resonator  and 
L  is  the  separation  between  the  two  resonators.  The  shape  of  the  electron  beam  is  contained  in 
the  sin  kya  factor,  so  that  pencil  beams,  annular  beams,  and  sheet  beams  can  be  considered.  The 
quantity  q  is  known  as  the  bunching  parameter  in  the  gyroklystron  literature,3  and  serves  as  a 
measure  of  the  phase  bunching  of  the  electrons  at  the  entrance  of  the  output  resonator.  As  will 
be  shown  later,  the  linear  efficiency  is  optimized  for  q  ~  2,  although  somewhat  higher  values  are 
required  in  the  nonlinear  regime. 
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B.  Nonuniform  Magnetic  Field 

In  the  existing  QOG  magnet,  the  field  is  nonuniform  for  much  of  the  region  proceeding  the 
output  resonator.  The  effect  of  this  nonuniformity  on  the  bunching  parameter  is  calculated  in 
this  section. 

If  the  magnetic  field  is  not  uniform  in  the  prebuiiching  resonator,  the  expression  for  Apj.  must 
be  modified.  The  non  relativistic  cyclotron  frequency  is  now  a  function  of  position 

fl(*)  =  flo  +  «l(*),  (12) 


where  6Sl(z)  is  a  small  quantity  compared  to  Qq  which  contains  the  nonuniformity  of  the  magnetic 
field.  The  change  in  perpendicular  momentum  is  now 


A. P±  =  sin  kyg  cos  [0O  +  +  Sil(z)  -  -)0u)z}dz.  (13) 

Z  J  -  oo  Pz\z)  Pz(z) 

It  is  assumed  that  the  energy  of  the  particles  is  constant  with  respect  to  position.  This  should  be 
a  good  approximation  for  the  purpose  of  calculating  the  momentum  change  in  the  first  resonator 
since  the  electric  field  is  small  and  the  change  in  particle  energy  is  small  (A7/70  <C  0.01). 

The  superconducting  magnet  currently  used  in  the  Naval  Research  Laboratory's  QOG  exper¬ 
iment  has  a  magnetic  field  with  a  positive  taper  for  much  of  the  region  preceeding  the  power 
resonator.  A  plot  of  the  magnetic  field  versus  position  is  shown  in  Figure  2.  The  slope  of  the 
magnetic  field  near  the  first  peak  is  approximately  0.9  kG/cm  at  the  position  z  =  - 10  cm. 
Th"  magnetic  field  may  be  written  B  =  Bo  +  6B(z),  where  6B(z)  -  0.92.  Here  2  is  given  in 
centimeters,  the  magnetic  field  is  in  units  of  kG,  and  z  =  0  is  referenced  to  the  center  of  the 
prebunching  resonator.  This  expression  may  be  substituted  into  Eq.  (13)  to  calculate  Apj..  The 
parallel  momentum  of  the  particle  is  also  a  function  of  axtal  position  aue  iu  the  variation  of  the 
magnetic  field.  Since  p\/B  is  an  adiabatic  invariant,  the  change  in  parallel  velocity  caused  by 


the  nonuniform  field  is 


6vt(z) 


-v2106B(z) 
2vzo  Bo 


(14) 


The  parallel  momentum  can  now  be  expressed 


P*(*)  =  Pzo[l  - 


<>26B(z) 
2  Bo  J 


(15) 


The  preceeding  expressions  for  E(z),  p?(z),  and  6Q(z)  can  now  be  substituted  into  Eq.  13.  Al¬ 
though  it  is  impossible  to  Derform  this  integration  analytically,  it  is  straightforward  to  implement 

numerically. 
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Typical  results  of  the  integration  of  Eq.  (13)  are  shown  in  Figure  3.  The  solid  curve  corre¬ 
sponds  to  the  change  in  perpendicular  momentum  for  a  particle  which  traverses  the  prebunching 
resonator  with  a  uniform  magnetic  field.  The  broken  curves  correspond  to  magnetic  field  gradients 
in  the  prebunching  resonator  of  0.5  and  0.9  kG/cm.  There  is  a  shift  in  gyropbase  angle  between 
the  uniform  and  nonuniform  cases  which  has  been  suppressed,  since  the  entrance  phase  into  the 
prebuncher  (#o)  is  arbitrary.  The  variation  in  gyrophase  shift  is  smooth  in  the  detuning  range 
between  ± 5 %.  All  of  the  plots  are  normalized  to  the  change  in  perpendicular  momentum  for  the 
uniform  field  for  zero  detuning. 

The  calculation  is  performed  for  a  frequency  of  85  GHz,  a  beam  voltage  of  70  kV,  and  o  =  1.5, 
where  a  =  vj_ / vjj .  The  mirror  separation  is  8  cm  and  the  radius  of  curvature  is  20  cm,  which 
correspond  to  the  values  of  several  available  mirrors.  It  is  seen  that  the  first  resonator  must 
be  driven  at  a  frequency  close  to  the  relativistic  cyclotron  frequency  for  the  prebunching  to  be 
effective.  As  the  magnetic  field  gradient  is  increased,  the  effectiveness  of  the  bundling  decreases. 
The  degradation  in  bunching  is  approximately  25%  for  a  magnetic  field  taper  of  0.9  kG/cm  for  the 
parameters  listed  above  at  zero  detuning.  The  radiation  waist  for  this  resonator  is  0.94  cm.  so  that 
the  magnetic  field  changes  approximately  ±5%  over  4  waist  radii.  The  effect  of  the  nonuniform 
field  can  be  made  smaller  by  decreasing  the  radiation  waist  in  the  prebunching  resonator. 

Figure  4  shows  the  results  of  the  calculation  of  A p±  for  a  =  1.0.  Although  the  shape  of 
the  curves  is  similar,  the  smaller  transverse  energy  results  in  a  smaller  change  in  perpendicular 
momentum.  Since  the  Apx  scales  as  1/p,,  the  ratio  between  the  two  curves  is  approximately  1.27 
tor  a  beam  voltage  of  70  kV.  The  bandwidth  assodated  with  the  prebunching  is  somewhat  larger 
for  lower  a.  This  information  is  important  because  it  is  difficult  to  get  an  accurate  estimate  of 
the  electron  pitch  angle  in  the  experiment.  A  spread  in  a  in  the  first  resonator  will  introduce  a 
spread  in  the  bunching  parameter,  therefore  the  nonlinear  effidency  should  be  calculated  for  the 
resultant  distribution  in  q. 

III.  Effect  of  Velocity  Spread  and  Nonuniform  Magnetic  Field  in  Drift  Region 
A.  Linear  Efficiency  in  Uniform  Magnetic  Field 

One  quantity  of  interest  is  the  small  signal  efficiency  when  a  uniform  magnetic  field  is  consid¬ 
ered.  The  energy  that  the  partide  gives  up  to  the  radiation  field  in  the  second  resonator  can  be 
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written 


AW  =  e  J  vxExdt  =  e  j ( px/pz)Exdz . 


The  expressions  for  the  electric  field  and  the  perpendicular  momentum  are  those  from  linear 


theory 


Px  -  P±  cos  =  Px  cos  (wt  +  0) 


4 

£r  - -Eo:  sin  fcys  exp  [ - 5-]  cos(ud  +  </>o).  (IS) 

to02 

Here,  </>o  is  the  phase  difference  between  the  rf  fields  in  the  two  resonators  and  the  subscript  2 
denotes  the  second  resonator.  Integrating  Eq.  (16)  and  averaging  over  the  gyroperiod  yields 

A„,/  a  A  ewoiVxEoipx  (fi0-70u ;)2wg2rn7,  .  , 

&W(yg,6 0)  =  - - - exp  [ - 73 - — — ]«n  kyg 

zPzO  4Pz0 

xcos[ — - — ^ - 2°^j,  L  -  4>o  +  0o  +  q  sin  kyg  cos  fluJ-  (19) 

PzO 

This  expression  for  A W  may  be  integrated  over  a  uniform  distribution  in  entrance  phase  #0  end 
an  arbitrary  distribution  in  the  y  direction  f(yg). 

<  AW  >  =  f-M'f  ^  ~ fa  -  ~  ^h}F{q)  ,20) 

^PzO  4p*0  pz0 


where 


E(g)  =  J  ^1(9  sin  Osin  £/v(0<^- 


The  function  F(?)  depends  upon  the  shape  of  the  beam,  where  f  is  the  variable  of  integration 
in  the  y-direction.  For  example,  a  pencil  beam  gives  F(q)  =  J\(q),  so  that  the  Unear  efficiency 
maximizes  for  q  —  1.84.  It  can  be  seen  that  the  linear  efficiency  of  the  gyrolystron  is  maximum 
for  zero  detuning 

— )-0.  (22) 

U>  U>lQ 

This  is  in  contrast  to  the  single-cavity  gyrotron,  where  the  output  frequency  is  always  greater  than 
the  relativistic  cyclotron  frequency.  However,  both  the  gyrotron  and  gyroklystron  reach  optimur 
efficiency  for  similar  positive  detuning  values.10 


B.  Velocity  Spread  and  Nonuniform  Magnetic  Field 

A  large  value  of  bunching  parameter  is  desired  so  that  the  gyroklystron  will  operate  at  peak 
efficiency.  This  can  be  accomplished  by  increasing  the  strength  of  the  prebunching  field  and/or 


increasing  the  separation  between  the  resonators.  A  large  electric  field  in  the  input  resonator 
will  give  a  large  perpendicular  momentum  kick,  which  will  introduce  a  sizeable  spread  in  energy. 
An  energy  spread  of  several  percent  is  undesireable  because  the  electrons  will  not  interact  effi¬ 
ciently  with  the  rf  fields  in  the  output  resonator.14  Increasing  the  separation  between  resonators 
makes  this  region  susceptible  to  velocity  spread  effects.  This  section  examines  how  the  gyrophase 
bunc  ing  is  affected  by  a  distribution  in  pitch  angle  and  a  nonuniform  magnetic  field. 

The  change  in  slow  phase  in  the  drift  region  is 

dz 

Aff  =  -  I  (u-  (23) 

Jo  T  vz 

If  the  magnetic  field  is  uniform  through  the  drift  region  and  the  electrons  are  monoenergetic,  the 
change  in  slow  phase  can  be  written 


fi  L 

AO  =  -(u - )— . 

7o»z 


(24) 


It  can  be  seen  that  a  spread  in  parallel  velocity  will  cause  a  distribution  in  gyrophases  for  electrons 
at  the  second  resonator  when  the  detuning  between  the  cyclotron  frequency  and  the  bunching 
frequency  is  not  zero.  When  this  spread  in  gyrophase  is  approximately  the  bunching  will 
completely  deteriorate  in  the  drift  region  and  the  device  will  operate  as  a  single- cavity  QOG. 

When  the  magnetic  field  in  the  drift  region  is  not  uniform,  the  cyclotron  frequency  and  parallel 
momentum  can  again  be  written  ft(z)  =  fio  +  6Sl(z)  and  pz(z)  =  pzo  -  dpz(z).  It  is  assumed  that 
7  =  7o  in  the  drift  region  for  the  present  discussion,  since  only  velocity  spread  of  the  beam 
electrons  is  considered  here.  The  change  in  slow  phase  in  the  drift  region  is  now 

fL  dz  ,  ,  fL  Sil(z)  dz  t  [L (  Q0^6vz(z) 

■4.  ' 


a  « =  -  /l  +  /‘  +  /V  - 

Jo  t>zo  7o  Jo  7o  Vzo  Jo  7o  vz0 


(25) 


7o  Vzo 

Here  the  second  order  term  proportional  to  6vz(z)6Sl(z)  has  been  dropped.  The  first  term  of  Eq. 
(25)  is  a  constant  with  respect  to  z  inside  the  integral.  The  other  two  terms  depend  upon  the 
nonuniform  field.  Define  the  average  variation  of  the  magnetic  field 

rL 


<6B  >=  ~  f  6B(z)dz. 
L  Jo 


(26) 


The  change  in  slow  phase  can  now  be  written 

Ad  = - (w - )  + - 

VzO  7o  rnioVzo 


<6B>  +Pr— (w  ~—)<6B> 
2tlo  7o 


(27) 
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The  technique  used  in  this  section  is  to  expand  A 8  in  terms  of  the  trajectory  pitch  angle  k, 


where  n  is  defined3 


vz  o  =  Vq  cos  K 


vio  =  Wo  sin  k. 


Now, 


.  .  L  ,fio  .  ,  eL  <6B>  L<6B>.Q, o  .  2 

AO  = - ( - w)  + - - - (—  —  u>)  tan2  k. 


vo  cos  k  '  70  '  m7o  tio  cos  k  2«o  cos  k  70 

The  change  in  gyrophase  angle  can  be  expanded  in  terms  of  6k  so  that 


(28) 

(29) 

(30) 


AO  =  Ai6k  +  ±A2(6k)2  +  A0o 

where  AOo  is  the  change  in  slow  phase  for  k  =  kq.  The  expansion  coefficients  are 

L  Q, q  eL 

A\  =  — (u> - )  sec  no  tan  kq - <6B>  sec/cotanKo 

wo  70  m70wo 

L  <6B>.  „  2  , 

+- - - — (w - )[2  -f  3  tan  KoJ  sec  kq  tan  kq 

Zvq  dq  70 


(31; 


(32) 


A2  =  —  (u  -  —  )sec«o[l  +  2  tan2  Ko) - <6B>  secK0fl  +  2  tan2  k0J 

wo  7o  m7o«o 

+  <^.>( u  -  — )  sec  Ko[2  +  13  tan2  ko  +  12  tan4  «o].  (33) 

-‘Wo  Do  70 

The  linear  efficiency  of  the  gyroklystron  is  a  convenient  measure  of  the  effect  of  velocity  sp  t  ad 
that  can  be  treated  analytically.  Thermal  spread  of  electron  velocities  will  principally  effect  e 
the  sine  factor  of  Eq.  (20)  because  its  argument  is  integrated  over  the  entire  drift  region.  Ti  > 
exponential  factor  is  due  to  the  force  bunching  in  the  output  resonator  where  the  integration  it 
performed  over  2tvo2.  Thus,  the  linearized  efficiency  with  velocity  spread  can  be  calculated  by 
averaging  over  the  sine  factor  if  the  drift  region  is  much  longer  than  the  radiation  waist  ( L  >  inoi , 
wo2)-  Assume  a  Gaussian  distribution  in  pitch  angle  k  =  «o  +  6k 

1 

/(8't)  =  7^“Pl-2?i-  (34) 

The  expression  for  <  AW  >  may  now  be  averaged  over  the  distribution  in  k. 

<  AW  >  =  ^etuo2\/*-Eo2 tan  ao  exp  [— (w  -  ~)2i^2/4r20]F(g)sin(<fo  -0) 

*  7o 

x(l  +  a4A2)_1/4exp[— ^A2a2/(1  +  aA  A\)\.  (35) 
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In  the  preceeding  equation, 


#  =  A0(*o)  +  ^  tan  1  Aid?  -  \-A\Aia4l(\  +  a4 A^).  (36) 

z  z 

The  small  signal  efficiency  of  the  gyroklystron  without  thermal  spread  is  obtained  by  setting 
a  =  0  in  the  proceeding  equations.  When  thermal  spread  is  considered,  the  coefficients  A]  and  A2 
modify  the  expression  for  efficiency  in  Eq.  (35).  These  coefficients  depend  upon  the  separation 
between  resonators  L,  the  mean  value  of  a,  and  the  magnetic  field  nonuniformity  <SB  >.  Define 
the  velocity  spread  coefficient 

C(a)  =  (1  +  a4Al)~1/4  exp  [— ~v4fa2/(l  +  a4A22)].  (37) 

This  coefficient  describes  the  degradation  of  the  linear  efficiency  due  to  velocity  spread  and  varies 
between  0  and  1.  The  result  of  the  calculation  is  shown  in  Figure  5  for  an  average  nonuniformity 
of  the  magnetic  field  of  4.1%,  which  is  the  value  for  the  QOG  magnet  for  L  ~  10cm.  The  mean 
value  of  pitch  angle  is  *0  =  1.0  ( a  =  1.56),  while  the  three  curves  correspond  to  pitch  angle 
spreads  a  =  0.025,  0.05,  and  0.1.  The  corresponding  values  of  A  a  are  approximately  5%,  10%, 
and  20%.  The  velocity  spread  coefficient  maximizes  at  a  frequency  detuning  of  approximately 
3.6%,  which  is  the  difference  in  cyclotron  frequencies  in  the  two  resonators. 

Figure  6  shows  the  results  of  the  calculation  of  the  thermal  coefficient  for  kq  =  0.8,  which 
corresponds  to  a  =  1.0.  The  full  width  at  half  maximum  for  the  three  curves  is  a  =0.025,  0.05, 
and  0.10,  yielding  similar  values  for  Act.  The  maximum  of  the  curves  is  shifted  slightly  to  the 
right  and  now  occurs  at  Au/u  =  3.7%.  The  point  where  the  effect  of  spread  is  minimized  depends 
only  weakly  upon  the  average  value  of  a,  so  that  fixing  the  frequency  detuning  doesn’t  constrain 
the  allowed  values  of  vx/u||.  The  main  dependence  of  the  optimum  frequency  detuning  is  on  the 
average  nonuniformity  of  the  magnetic  field.  The  magnetic  field  shape  can  be  varied  somewhat  in 
the  experiment,  which  allows  for  a  technique  to  tune  the  optimum  frequency  difference  between  the 
input  and  output  resonators.  Another  feature  of  Figure  6  is  that  the  bandwidth  of  the  interaction 
increases  with  decreasing  oo  for  large  veloaty  spreads.  Thus,  it  may  be  desireable  to  operate  in 
thi8  regime  to  lower  the  effect  of  thermal  spread. 

IV.  Estimation  of  the  Bunching  Parameter 

The  condition  for  starting  oscillation  in  the  first  resonator  with  an  annular  electron  beam  with 
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optimum  frequency  detuning  can  be  expressed2,4 

(38) 

In  the  above  equation  Pi,  is  the  input  beam  power  in  Watts,  d  is  the  mirror  separation,  rj,  is  the 
beam  radius,  and  0  is  the  ratio  of  the  particle  velocity  to  the  speed  of  light.  Considering  the 
annular  beam  is  important  because  some  electrons  pass  through  nulls  of  the  standing  fields  in 
the  resonator,  which  decreases  the  efficiency  of  the  interaction.  If  the  guiding  center  radius  rf,  is 
somewhat  greater  than  the  wavelength,  then  Jo(kri,)  ~  0  and  Eq.  (38)  can  be  written 

KQ  >  4.6  x  109(— )7o(7o  -  l)^2-.  (39) 

r01  P±  0 

As  an  example,  consider  a  resonator  with  a  mirror  separation  of  8  cm  and  mirror  radius  of 
curvature  20  cm.  For  a  frequency  of  85  GHz,  the  radiation  waist  is  0.94  cm.  An  electron  beam 
voltage  of  70  kV  corresponds  to  a  relativistic  mass  factor  70  =  1.137,  and  assuming  a  =  1.5  yields 
Pxo  =  0.40  and  0zo  =  0.266.  If  the  quality  factor  (Q)  of  the  resonator  is  1000,  then  the  threshhold 
current  of  the  prebunching  resonator  is  9.6  A.  The  Q  of  the  input  resonator  is  determined  by  the 
diameter  of  the  mirrors  and  the  size  of  the  coupling  hole. 

The  electric  field  in  the  bunching  resonator  determines  the  magnitude  of  the  bunching  param¬ 
eter.  A  fundamental  relation  between  the  energy  stored  in  the  resonator  and  power  dissipated  by 


losses  is 


Q  =  uWem/Pl, 


where  W„n  is  the  stored  electromagnetic  energy  and  Pi  is  the  average  power  lost.  For  a  Fabry- 
Perot  resonator,  the  stored  energy  can  be  written15 

=  £co  w^dEl  (41) 

Here,  €0  is  the  permeability  of  free  space  and  all  quantities  are  in  MRS  units.  For  the  resonator 
described  above  and  an  input  power  of  375  Watts,  the  electric  field  at  the  center  of  the  resonator 
is  £bx  =  1-67  x  10®  V/m.  The  total  power  available  from  the  EIO  is  1500  Watts,  so  the  value  375 
Watts  should  be  a  conservative  estimate  of  the  power  coupled  into  the  TEMoo  mode. 

If  the  drive  frequency  is  sufficiently  close  to  the  relativistic  cyclotron  frequency  in  the  first 
resonator,  then  the  exponential  term  in  Eq.  (11)  is  approximately  one.  When  the  effect  of  the 
magnetic  field  taper  is  included,  the  bunching  parameter  can  be  written 

q  ~  (0.75)v/5fwpxeEoitn0iL/2p^0c2. 
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(42) 


The  separation  L  between  the  two  resonators  is  10  cm  for  a  detuning  of  3.6%,  yielding  a  bunching 
parameter  q  =  1.8.  Tran  et  al.  have  shown  that  the  optimum  bunching  parameter  for  peak 
efficiency  with  a  pencil  beam  placed  on  a  field  maximum  is  approximately  3.10  The  optimum 
value  of  q  for  an  annular  beam  is  somewhat  greater  due  to  variation  of  electric  field  across  the 
beam.  Thus,  the  maximuir  'nching  in  this  example  is  somewhat  less  than  optimum. 

It  can  be  seen  from  Eq.  (11)  that  the  bunching  parameter  scales  as  Px/p],  so  that  reducing 
vs/vji  to  1.0  reduces  the  bunching  parameter  by  approximately  50%  for  a  beam  voltage  of  70  kV. 
This  operating  regime  (low  a)  is  of  interest  because  velocity  spread  effects  will  be  smaller  and  the 
gyroklystron  can  operate  at  higher  current  before  the  prebunching  resonator  oscillates 

V.  Conclusions 

The  slow  time  scale  equations  of  motion  have  been  linearized  to  study  a  quasioptical  gy¬ 
roklystron  with  a  nonuniform  magnetic  field  in  the  prebunching  resonator  and  drift  region.  The 
prebunching  resonator  must  be  driven  close  to  relativistic  cyclotron  frequency  for  the  prebunching 
to  be  effective.  A  tapered  field  in  the  input  resonator  is  shown  to  reduce  the  bunching  parameter 
by  approximately  25%  for  a  taper  of  ±5%  over  4  waist  radii.  This  effect  can  be  reduced  by 
decreasing  the  radiation  waist  in  the  prebunching  resonator. 

The  equations  of  motion  are  integrated  through  the  drift  region  to  examine  the  effect  of  velocity 
spread  coupled  with  the  nonuniform  magnetic  field.  It  is  assumed  that  the  drift  region  is  much 
longer  than  the  radiation  waist  in  the  two  resonators  so  that  the  predominant  effect  of  velocity 
spread  is  phase  mixing  of  the  beam.  For  the  experimental  parameters  chosen  here,  a  frequency 
detuning  of  3.6%  minimizes  the  effect  of  the  velocity  spread  for  the  magnetic  field  variation  of 
4.1%.  Thus,  the  difference  between  the  cyclotron  frequencies  in  the  two  resonators  should  be  fixed 
at  3.6%,  which  should  result  in  good  nonlinear  efficiency.  For  the  available  millimeter- wave  source, 
a  maximum  bunching  parameter  of  q  ~  1.8  is  obtained.  Varying  the  bunching  from  small  to  large 
values  should  allow  gyroklystron  operation  in  several  different  regimes  so  that  mode  selection  and 
efficiency  enhancement  can  be  investigated. 
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Schematic  diagram  of  a  two-resonator  quasioptical  gyroklystron 
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Fig.  5  —  Velocity  spread  coefficient,  defined  in  F.q.  (37),  vs  frequency  detuning.  The 
pitch  angle  distribution  is  Gaussian  with  mean  <i0  =  1.0  and  standard  deviation  a  =  0.025,0.05, 
and  0.1  for  the  solid,  dashed,  and  dotted  curves,  respectively. 
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